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Abstract 

We explore the possible occurrence of a-uo condensation in the 
Quantum Hadro-Dynamics (QHD), namely the Serot and Walecka 
model, finding that at the mean field level it corresponds to a critical 
value of the coupling constant g a = 8.828 and density kp = 207.2 
MeV/c, significantly below the standard value of QHD. 

PACS: 21.65.+f, 24.10.Cn, 24.10.Jv 

1 Introduction 

The Quantum Hadro-Dynamics (QHD) is a largely used model to describe 
nuclear systems within a coherent and covariant frame. Its parameters (cr 
and u masses and coupling constants) are tuned to reproduce the energy 
and density of the nuclear matter at the mean field level, where exotic or 
non-trivial phenomena are intrinsically forbidden. 

In going beyond the mean field, however, one could question above the 
stability of the nuclear matter against the occurrence of a a-u condensate. 

We wish to establish in this paper the limits the stability imposes on the 
model parameters. 
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2 General formalism 

To shortly remember the QHD scheme PQ, we introduce the lagrangian den- 
sity 

C = C N + C rT + C tJ + C I (1) 



with 



C N = ft-m)^ (2) 

= \{^f-\m 2 y (3) 

Cu, = —F^ + ^mlu^ (4) 

£i = 9aipcrtfj - gj^rfp^uf - ^a 4 a 4 (5) 



and 



F» v = d»uf - ffV . (6) 

We get rid of the non- vanishing a and u° fields in the vacuum by means 
of the shifts a = a' — a, — > u^' — cD^, the mean fields being determined 
by solving the classical equations of motion: 



mla + |yfT 3 = g a < ^ > (7a) 

™l& = 9u> < i^loip > (7b) 
u i = (7c) 



Their solutions read 



and 



fhP . 
m 2 , 



<P° (8) 



where p is the usual nuclear density 



3tt 

and 



R a = V 3\falg<rPa + V 8m a + 9a 4 ^p2 (H) 



with ^ 

P<y = -^r J ^yS H {k) . (12) 

In the latter definition Sh denotes a fermion propagator analogous to the 
nucleon propagator in the medium 

s^k) = u ' + m { - kp) + e ^ F - k) 1 

2E k \k -E k + ir] k - E k - ir] k + E k - %r\ 

but with the nucleon mass replaced by an effective one, defined by 

m* = m — a . (14) 

Eqs. (jnj) and (fTT|) need to be solved self-consistently in a. Its knowledge 
immediately provides the effective nucleon mass. It is interesting to note 
that a self-consistent calculation realizes the prediction of Lee and Wick j2] 
that for large kp the nucleon mass is vanishing. 




3 The Random Phase Approximation (RPA) 

It has recently been observed [3] that the presence of a nuclear medium couples 
the a to the th component of the u, so that the a and u propagation in the 
medium will be described by two coupled Dyson equations. We can write 
them (in momentum space) as a single equation in a 5-dimensional space in 
the form 

D (5) = £,(5) + jD (5) n *(5) jD (5) ( 15 ) 

where and n*( 5 ) are 5x5 matrices having the structure 

( 5 ) _ fD 0a o \ r) (b)_(E>*a DL\ n *(5) _ /n* CT n*^\ . , 
" V o dz) u - \Di a dz) ~ \nz nrj ' [ b) 

The free meson propagators are defined as 

Doa(q) = 9*— 2 \ , . g a , (17) 



DZ(q) = -qA g u = -(.</"" -^) A„.(v). (is) 



3 



The polarisation propagators are constrained by current conservation. Thus 
they must take the form 



KM) = KM)w = u v ( q )w , 



with 
and 

with 
and 



1, 



q°q l 



K V M) 





90 % Yi L 
JqP" 






WJ 



iT 
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q l q J 



nl(g) . 



The quantity T\* aa (q) has no tensor structure. We shall define 

n s = n* 



(19) 
(20) 

(21) 

(22) 
(23) 

(24) 



in order to harmonise the notations. 

Further, since q^D^ = q^D^ = and q^D^ = 0, as Dyson equation 
implies, we must also have 



and 



D 


v = 

era; 


LUU 




D L 










w 


DL w + ^ T 0- 


MV 



(25) 



(26) 



The polarisation propagators needs to be approximated in some way. It 
has been proved in [3] that at the level of mean field in a bosonic space they 
must be replaced with their zero-order approximation. Thus we replace IT* 
with IT , with components (traces refer both to spin and isospin) 



f d 4 p 

IW(g) = -iTr J -j—^S H {p)S H {p + q) , 



(27) 
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and 



r d 4 p 

n o^(?) = -iTr J j—^S H (p)S H {p + 



(28) 



+ q)Y ■ 



(29) 



The th order propagators Il have been extensively studied in ref. 
jS], where their analytical representation is given. Here we observe that they 
can be written as 



n, 



4(4m 2 - q 2 )U° + 8T 
16m£3 y 

I I 2 

16£3 L -4|q| 2 n° + 8^X 
16£3 T + 8g 2 n° - 161 , 



(30) 
(31) 

(32) 

(33) 



where 

n° = 
& v = 



d 3 p 1 6(k F — p) 
(2^2^(96 + ^-^ 

d 3 p t° 6(k F - p) 
(^^(gb + ^-E^ 

d 3 p (ff 9{k F -p) 
(2vr)3 2£ p (g + £ P ) 2 - £ 2 +q 



+ (90 



-go) 



Po=-E P 



+ (9b < — * -Qo) 



Po=E P 



+ (9o 



-Qo) 



(34) 
(35) 
(36) 



Po=E p 

d 3 p |p| 2 |q| 2 — {p ■ q) 2 9(k F — p) 
^q] 1 (90 + £ p ) 2 - 



p+q 



Po=E p 



d 3 p 6{k F — p) 
(2tt)3 2E p ■ 

In the above we have introduced the transverse vector 



f = - 



p_q 

q 2 



+ (go < — ► -90) , 

(37) 
(38) 



(39) 



while 

% = <?o + iv sign(g ) (40) 

accounts for the right analytical determination near the cuts. 
We can now solve the Dyson equation ()15j) . finding 

D aa = D a 1 (41) 

i - uXd,Md„ 



where 



1 2 



*0 ±J u) 1L 

2 



D L =W\\ D u -4 ~ (42) 

ii 2 2 ~ 

D V = _]S^ DmU V Du = _y DaU V D L (43) 
Q |q| 

D T = \ Dou T , (44) 



^ = ; D °°s (45) 



1 - D 0(r Ul 



2 



|q| 2 



Observe that the transverse part of the cu-meson is completely decoupled 
from the longitudinal propagation, that instead involves a and u°. 



3.1 The (j-oj condensate 

We consider first a simplified model where only the a meson exists. Since 
the cr-exchange is attractive, a a condensation is expected under certain 
conditions. In this model the RPA a propagator (|41j) reads 

and a a condensate arises if the denominator of the above equation vanishes 
at a certain |q| and at q = 0, i.e., if the equation 

1 - £W(|qUo = 0)n*(|qUo = 0) = (48) 
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admits some solution in |q|. 

Let us now study the function £) 0o .(|q|, q = 0)IIo (|q| , g — 0). From 
eq. (fTTj) we get 



-DoadqUo = 0) 



9l 



|q| 

while an easy calculation, using the explicit expressions given in [5], provides 



n 5 (|qUo = o) 



4m 2 + |q| 



2vr 2 |q| 



|q| log 



k F + E F 



rn 



+ E F log 



2k F + |q| 



2k t 



-\l Am 2 + |q| log 



E F \q\ + k F \ /Am 2 + |q 



i?i?|q| — k F \ I Am 2 + |q 



(49) 



A simple check shows that this function is regular at |q| = 2k F , and from 
the above it immediately follows that 



lim D 0a (\q\,q = 0)n o 5 (|qU = 0) = ^_^\ og ^±^ > 



|qH0 



m 2 it 2 



(50) 



On the other hand for large |q|'s LTq (|q|, qo — 0) has a finite limit so that 



D 0a (\q\,q o = 0)Ul(\q\,q = 0) 
As a consequence eq. 



|q| 



|q|— »oo 

has certainly a solution provided 



gl 2m 2 k F + E F 

log > 1 



m 2 7r 2 



rn 



(51) 



(52) 



The physics contained in this conclusion was expected: a sufficiently large 
attraction produces a a condensate. Numerically at the normal nuclear den- 
sity [kp = l.SGfur 1 ) and with m a = 550 MeV used in QHD [T] we get for 
g a the critical value g a = 2.47, well below the value of the QHD, namely 
^ = 9.573. 

Actually the w-meson exchange, being repulsive, could prevent such an 
occurrence. A simple algebra shows that D aa , D L and D v , as given in eqs. 
(I41|) . (|42j) and (J43|) . have the same denominator, namely 

S)l = [1 - D 0a U s ] [1 - (C - 2)D 0uj U L ] - (C - 2)D 0a D 0uJ (U v ) 2 (53) 
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(the index L reminds us that we are concerned with the longitudinal propa- 
gation), where we have introduced the shortcut 



c = 



|q| 


2 + <i 




q 


2 



C-2 = 



|2 ' 



(54) 



The occurrence of a phase transition is now signalled by the existence of 
solutions of the equation 

®l = . (55) 



Proceeding as before we find 



n y (|q|,g = o) 
n L (|q|,g = o) 



m 



(27r)2|q| 
1 

'67r 2 |q| 



4A: F |q| + (Ak 2 F - |q| 2 ) log 



2k F + |q| 



2\q\k F E F + 2|q|(3m^ - |qf ) log 



2k F — |q| 

k F + E F 



m 



(56) 
(57) 



+ E F (AE F -3\qf)\og 



2k F + 


|q| 


2k F - 


|q| 



yj 4m 2 + |q| 2 (2m 2 — |q| 2 ) log 



E F \q\ + kp^Am 2 + |q| 



E F \q\ — k F \J Am 2 + |q| 



together with the relevant limits for |q| — > 

2mk F 



lim n K (|q|, = 0) = 



|q|-0 



7T Z 



lim n L (|q|, g = 0) = <^ kpEp + m 2 log 



|q|->0 



7T Z 



A; F + E F 



m 



(58) 
(59) 



both limits being negative, and for |q| — > oo 

1 



n^(|q|,g = o) 

iqi 

n L (|q|,go = 0) — > constant 
Moreover in the limit ^ ^ we find 

C = i, C-2 = -i 



(60) 
(61) 
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Thus at q = 

£l | 9o=0 = [1 - D 0a U s ] [1 + D 0u U L ] + D 0a D 0u) {U v ) 2 . (62) 
which implies, using (flgj) . (JoTj) . (jHOJ) and (pTTj) 

lim S) L (|qUo = 0) = l . (63) 

|q|-»oo 

Thus a condensed state will occur if 



lim S L (|q|,go = 0) <0 . (64) 

|q|->0 

Let now consider the above limit as a function of k F : 

<P(k F ) = lim D L (|q|,g = 0) . (65) 

|q|->0 

From (J5T]) . (p)H|) and (JSTJj) we derive, for low kp, 

<j>{k F ) = 1 - ^ f 4 - 4 V + , (66) 

while at large k F we find the asymptotic behaviour 

«*,) * -^44» 2 4 lo« — + . (67) 

Hence <p{k F ) goes to 1 for low kp and to — oo for large kp and thus a critical 
value for /cp surely exists. The problem of how and when this critical point 
is reached and this can be addressed only numerically 

Before exploiting the calculations, we observe that up to now we have 
not accounted for self-consistency in determining the nucleon mass, accord- 
ing to eq. ()14j) . The nucleon effective mass m* is displayed in fig. ^ and 
compared with a non-self-consistent calculation (i.e., with Sh replaced by So 
in eq. (fT2jO. The plot shows, as pointed out by Lee and Wick many years 
ago|2], that the self-consistent effective mass tends to vanish at large den- 
sities. Thus the large kp limit (J67|) must be evaluated at m = 0, namely 



<j>{k F ) ~ 44 > • (68) 
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Figure V. Nucleoli effective mass. Solid line: full self-consistent calculation; 
dashed line: self-consistent calculation with = 0; dotted line: non- self- 
consistent calculation, dash-dotted line: non- self-consistent calculation with 
a 4 = 0. 

On the other hand, at low hp the effective mass tends to the bare one 
and with the standard parameters of QHD, namely g a = 9.573, m a = 550 
MeV, g w = 11.67 and = 783 MeV, 4>{kp) start decreasing if 



Thus one may reasonably argue that the function may have a minimum below 
0. An explicit calculation confirms this guess, as shown in fig. EJ 

Actually, to be reasonably sure that no condensation arises, we could 
require 

L=o> 

that corresponds to a limiting value g a = 8.197, to be compared with g a = 
9.573 as given in [Tj. 

A more detailed search for limiting conditions requires a numerical cal- 
culation. In fig. 01 we show the behaviour of (ft for five different values of g a 



m 2 , 
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Figure 2: Plot of the function 4>{kp) with a self-consistent effective mass 
(solid line) and the bare mass (dashed line). The last is evaluated with the 
standard values of QHD. 



ranging from g a = 10 (lowest solid line) to g a = 8 (upper solid line). It turns 
out that the critical value for the existence of a a- condensate is 

9<TcAt = 8.828 . (69) 

It is clear nevertheless that close to the critical value precursor phenomena 
(like in pion condensation) are to be expected, thus a lower value for g a 
should be preferred. 

The above states that at g a = 8.828 there is (as the curves in fig. El 
suggest) only one value of kp at which the cr-condensation occurs. Numeri- 
cally this happens at k c p lt = 207.8 MeV/c. Above this value we see that two 
solutions of the equation 4>{kp) = exist. Thus below the lowest solution we 
find no cr-condensation, and the same occurs above the higher one. In fig. 0] 
we have plotted the function Dl| 90= o for three different values of hp. This 
shows how below the lower value and above the higher one no phase tran- 
sition seems to occur. Actually this statement is only formal: indeed when 
a phase transition is found at the lower critical value of hp, then what hap- 
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Figure 3: The function </>(&f) for different values of g a , ranging from 10 
(lowest solid line) to 8 (upper one) with a spacing of 0.5, with self-consistent 
effective mass. The dashed line corresponds to a bare nucleon mass. 



pens for higher kp becomes unpredicable in the present formalism. On the 
other hand since no evidence exists about this phase transition, the existence 
of another critical value of kp is simply of mathematical, but not physical, 
interest . 

4 Conclusions 

The QHD was taylored to describe the static properties of the nuclear matter 
at the mean field level, i.e., at the th order in the loop expansion. Within 
this approximation no phase transition can arise. 

We can go beyond, however, and consider an expansion in bosonic loops 
only (boson loop expansion, in short BLE), as described in |3J Ej. There 
it was proved that the RPA series is just the mean field of the BLE: it 
contains in fact fermionic loops but not bosonic ones and in principle it may 
originate phase transitions, as is well known from the old discussions about 
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Figure 4: The function DL\ qo =o for different values of k F . Solid line: 
k F = 127.8 MeV/c, dashed line: k F = 207.8 MeV/c, dotted line k F = 287.8 
MeV/c. 



pion condensation. 

In this paper we have shown in fact that at BLE-mean field level the 
Serot and Walecka model would predict the occurrence of a g-uj condensed 
phase. This outcome however does not destroy the whole apparatus of QHD: 
it simply states that the mean field level is not able to provide a stable ground 
state at least above some critical values of k F . It is remarkable, nevertheless, 
that these critical values lie below the normal nuclear density. 

To overcome this instability we need to perform a higher order calculation 
of the binding energy (at the order of two bosonic loops). Since the value 
of the parameters in this model have been fixed at the mean field level, in 
going beyond we need a reparametrization, that hopefully should lead to a 
non-critical choice. 

To conclude we observe that in the present model we have studied the 
possible occurrence of a condensed state as g a varies. We could as well modify 
the a mass, as the relevant parameter is their ratio. Actually this is ininfluent 
as far as we look at the existence of a critical point, however, it matters when 
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we consider the momentum q where the condensation occurs. 
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